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^ ' Abstract 

We discuss single transverse spin asymmetries in the Drell-Yan process originating from 
so-called gluonic poles in twist-three hadronic matrix elements, as first considered by Qiu 
and Sterman. Even though time-reversal invariance is not broken, the effects of such poles 
cannot be distinguished from those of time-reversal odd distribution functions. We show 
the connection between gluonic poles and large distance gluon fields, in particular we focus 
on boundary conditions. We identify the possible single spin asymmetries in the Drell-Yan 
process. 



1 Introduction 



In the standard description of the Drell-Yan process (DY) in terms of distribution functions time- 
reversal symmetry implies the absence of single spin asymmetries at tree level, even including 
order 1/Q corrections [l]]. Additional time- reversal odd (T-odd) distribution functions (DFs) 
are present when the incoming hadrons cannot be treated as plane-wave states. This may 
occur due to some factorization breaking mechanism ||. We will show that, even apart from 
such mechanisms, the contributions of T-odd DFs may effectively arise due to the presence of 
so-called gluonic poles attributed to large distance gluon fields. The gluonic poles appearing 
in the twist-three hadronic matrix elements ^, |5[] together with imaginary phases of hard 
subprocesses effectively give rise to the same single spin asymmetries as T-odd DFs, but without 
a violation of time-reversal invariance. This is the origin of the single spin asymmetry of Ref. 
|J. For a detailed account on these matters see J7j. 

1 Invited talk presented at the workshop 'Deep inelastic scattering off polarized targets: theory meets experi- 
ment', DESY-Zeuthen, September 1-5, 1997 
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2 The DY process in terms of distribution functions 



We employ methods originating from Refs. || |9|, [K| |TT|, |12|, [13|, [14|] in order to describe the 
soft (non-perturbative) parts of the scattering process in terms of correlation functions, which are 
(Fourier transforms of) hadronic matrix elements of non-local operators. We restrict ourselves to 
tree-level, but include 1/Q power corrections. The asymmetries under investigation are loosely 
referred to as 'twist-three' asymmetries, since they are suppressed by a factor of 1/Q, where the 
photon momentum q sets the scale Q, such that Q 2 = q 2 . We do not take Z bosons into account, 
since the asymmetries are likely to be negligible at or above the Z threshold. 

The Drell-Yan process consists of two soft parts (depicted in Fig. [I] for the leading order) and 
one of them is described (up to order 1/Q) by the quark correlation functions $ and &a, and the 
other soft part by the antiquark correlation functions, denoted by $ and The quark-quark 




Figure 1: The leading order contribution to the Drell-Yan process 



correlation function, 

QiiiPuSxp) = J ^e^(P 1 ,S l \^ j m i (z)\P 1 ,S 1 ), (1) 

is a function of the momentum and spin vectors Pi, Si of the incoming hadron (spin-1/2), 
with Pi • Si = 0, and the quark momentum p. The hadron momentum Pi is chosen to be 
predominantly along a light-like direction given by the vector n + . Another light-like direction 
n_ is chosen such that n + -n_ = 1; both vectors are dimensionless. The second hadron is chosen 
to be predominantly in the n_ direction, such that Pi • Pi — 0{Q 2 ). We write p ± = p ■ n T and 
approximate the parton momentum p « xPi+px and the polarization vector Si ~ \iP\/Mi+S 1T . 

We will consider the case where one integrates over the transverse momentum gr of the 
photon. One then only encounters correlation functions integrated over all but the leading 
component, such that they are functions of the light-cone momentum fractions (e.g. x) only. So 
we consider the partly integrated quark correlation functions 

^(x) = J — e^(Pi,5i|^.(0)^(An_)|Pi,5i), (2) 

n^,y) = [ ^p-e^j^iPi, 5 1 |^.(0)^(7 7 n_)^(An_)|Pi, Si). (3) 
J J In Zn J 
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The anti-quark correlation function is defined as 



Hx) = I ^-e-^(P 2 ,S 2 m\n + )m\P2,S 2 ) 



(4) 



and the function & A (x, y) is defined analogously. Moreover, the inclusion of path-ordered expo- 
nentials, such as, 

£(0, An_) = V exp (-ig dz^ A^z)] , (5) 

which are needed in order to render the correlation functions gauge invariant, is implicit. 

In the expression of the hadron tensor the fermion propagators appearing in the hard part of 
the subleading contributions (cf. Fig. |2|) is approximated like (neglecting contributions that will 
appear suppressed by l/Q 2 ) 

j> x - j + m _ j + x-y 
(Pi _ o) 2 - m 2 + ie Q\/2 x - y + ie 

Hence, a zero-momentum gluon (x = y) is always accompanied by an on-shell quark propagator, 




4 — fW^ VIP 



Figure 2: Subleading order contribution to the Drell-Yan process 



where we note the following: 



dyn&y) - y ■ / dyn&v)- 



, ., , . (7) 

x — y + ie J 

For the correlation functions $ and ^ we need up to order l/Q the following parametriza- 



tions in terms of distribution functions [14| : 



$(x) = -[/i(x)A + g 1 (x)\ 1 >y 5 p 1 +h 1 (x)>y 5 p iT Pi 



+ 



Mi 

~Y 

Mi 



e(x)l + g T (xh5fiiT + h L (x) y7s [fi+, i-\ 



G A (x, y) idfSvrpft + G A (x, y) S^y s JPi 



(8) 



+ H A (x,y)\ ll5l %f 1 +E A (x,y) 1 2P 1 



(9) 
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where €j? 



e a ^ v n +a n-p. 



The parametrization of is consistent with requirements following from hermiticity, parity 
and time-reversal invariance, 



&(Pi,Si,P) =7o$(i 3 i,£i;p) 7 o 
$(P 1 ,5 1 ;p)=7o$(A,-5 1 ;p) 7 o 



[Hermiticity] 
[Parity] 

[Time reversal] 



(10) 

(11) 
(12) 



where p = (p°, —p), etc. For the one-argument functions in Eq. @ it follows from hermiticity 
that they are real. Note that for the validity of Eq. ( p"2"D it is essential that the incoming hadron 
is a plane wave state. For Q A hermiticity, parity and time-reversal invariance yield the following 
relations: 

[^ a A (P l ,S 1 ;p ll p 2 )] j = 7o^(P 1 ,£' 1 ;p 2 ,p 1 )7o [Hermiticity] (13) 
n{Pi,Si;Pi,p 2 ) =7o$Aa(A,-£i;Pi,p 2 )7o [Parity] (14) 

[^{P 1 ,S l ;p l ,p 2 )}* = l5 C$ Aa (P 1 ,S 1 ;p 1 ,p 2 )C^ 5 [Time reversal] (15) 

Hermiticity then gives for the two-argument functions in Eq. @ the following constraints: 



G A (x,y) = -G* A (y,x) 
E A (x,y) = -E* A (y,x) 



G A (x,y) = G A (y,x), 
H A (x,y) = H A (y,x). 



(16) 
(17) 



Hence, the real and imaginary parts of these two-argument functions have definite symmetry 
properties under the interchange of the two arguments. If we would impose time-reversal invari- 
ance all four functions must be real and G A and H A are then symmetric and G A and E A are 
antisymmetric under interchange of the two arguments, such that at x = y only G A and H A 
survive. 

In the remainder of this section we do not impose time-reversal invariance and hence allow 
for imaginary parts of these functions. In addition, the following (T-odd) one-argument DFs 
then appear: 



T-odd 



f T (x)£T 'SiTpjTv ~ e L (x)\ii-i b + h(x)- [fl + , flJ\ 



The two-argument functions and the one-argument functions are related by the classical 
e.o.m. ((ift — w)ip = 0), which hold inside hadronic matrix elements Using a similar 

parametrization for $J (defined like <&% but with gA9p replaced by iDj) as in Eq. @, one has 
the following relations O, H: 



dy ReG D (x,y) + ReG D (x,y) 
dy lmGi)(x,y) + lmGi)(x,y) 
dy 2 Re H D (x,y) 
dy 21m H D (x,y) 
dy 2ReE D (x,y) 
dy 2lmE D (x,y) —ixh(x). 



m 



2xg T (x) - 2—hi(x), 
2ixf T (x), 



in 



xh L (x) - j^9i(x), 
—ixe^x), 



in 



(19) 
(20) 
(21) 
(22) 
(23) 
(24) 



4 



From this (and iD a = id a + gA a ) we see that the (T-odd) imaginary parts of the two-argument 
functions are related to the T-odd one-argument functions, as one expects. So if time-reversal 
invariance is imposed, the imaginary parts of the e.o.m. Eqs. (|20|), ([22|) and ([24]) become three 
trivial equalities. We like to point out that if one integrates Eqs. ( Pf ) and fl2D| ) over x, weighted 



with some test-function cr(x), one arrives at the sum rules discussed in [I3L [15 



3 Gluonic poles and time-reversal odd behavior 

We are interested in the behavior of the quark-gluon correlation function $^ in case x = y, when 
the gluon has zero-momentum. For this purpose, we define (a is a transverse index) 

$%(x,y) = J ^^e a V^(P,5|^.(0)F +a (^_)^(An_)|P,5) (25) 

and F pa (z) = ~ [D p (z), D a (z)}. This matrix element has the same hermiticity, but the opposite 
time-reversal behavior as 

[<f> F {x, y)]* = -7sC $ Fa {x, y) CS [Time reversal] (26) 

and we will parametrize it identically with help of functions called G F (x,y),G F (x,y), H F (x,y) 
and Ep(x,y), noting that time-reversal implies that G F and Ep are symmetric and thus may 
survive at x = y (in contrast to Ga{x, x) and Ea{x,x)). In the gauge A + = one has F +a = 
d + A T and one finds by partial integration 

(x-y)n{x,y) = -m(x,y). (27) 

If a specific Dirac projection of $> F (x, x) is nonvanishing, then the corresponding projection 
of ^%(x, x) has a pole, hence the name gluonic pole. An example is the function T(x,St) = 
yrTr \® F (x, x) e T0a S^flJ / P + = 2mMS T G F (x,x) discussed by Qiu and Sterman in Ref. ||. 

In order to define Eq. at the pole, one needs a prescription, which is related to the choice 
of boundary conditions on Ap(r] = ±00) inside matrix elements. Possible inversions of F +a = 
d + Ap are: 



/oo 
dz~ 6{z~ -T 1 n_)F +a (z-) 
-00 

/oo 
dz~ 9(r]n^ - z~) F +a (z 
-00 



Ap{oo) + A%{-oo) 1 /-^ 



2 2 
One can use the representations for the 9 and e functions 



dz~ e(z~ -rjn_)F +a {z-). (21 



±a ^= ^=/i p V- (29) 



k =F it 



to obtain 

A(oo) 



n(x,y) = S(x-y)^ A(oo) (x)+ \ $ F (x,y) (30) 
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S(x-y) $3(-oo)0>0 + 



x — y — it 



5(x - y) 



A(oo) 



f x) + $ 



A(-oo) 



X 



+ P- 



—% 



x-y 



(31) 



where 



5(x-y) $S(± 



oo)ij 



^e a V^-*)(P, 51^(0)^(77 = ±oo)^(An_)|P, 5). (32) 

Z7T Z7T J 



So Eq. (|3T|) shows the importance of boundary conditions in the inversion of Eq. (|27|), if matrix 
elements containing Ay (77 = ±00) do not vanish. When such matrix elements vanish (implicitly 
assumed in |j]]) the pole prescription does not matter. Also one obtains 



2vr $%(x,x) 



(33) 



which shows the relation between the zero-momentum quark-gluon correlation function and the 
boundary conditions. 



The behavior of 3>5l(±oo) ( x ) un der time-reversal is: 

*A(±oo)(^) = 75C $A( T oo)a(x) C^ 5 . 



(34) 



This relation implies that time-reversal invariance only allows for symmetric or antisymmetric 
boundary conditions. Both situations (if nonvanishing) lead to a singularity in ^%(x, y) at 
the point x = y, but only the antisymmetric case will be called a gluonic pole. The delta- 
function singularity in the case of nonvanishing symmetric boundary conditions will contribute 
to the functions Ga(%, x) and Ha(x,x) and hence, to T-even DFs. This would only affect the 
magnitude of double spin asymmetries. This case is also less interesting, because $>p(x, x) = 0. 

We like to point out that so-called fermionic poles play a role in off-forward scattering, such 
as prompt photon production [[U], Here a gluonic pole gives rise to an asymmetry 

proportional to T(x, ST)g{x) (see Fig. |3|). Fermionic poles do not contribute in case of DY to 
this order. 




Figure 3: A diagram yielding a single transverse spin asymmetry in prompt photon production 
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4 Effective T-odd distribution functions 



To study the effect of gluonic poles we consider the (nonvanishing) antisymmetric boundary 
condition &%oo){ x ) = ~ ®A(-oo)( x )i which implies 

= *% 00) (*) > (35) 
n\±oo)(x) = -l 5 C$ A{±oo)a (x) Ct T5 . (36) 

In the calculation of the cross-section one always encounters the pole of the matrix element 
(in this case in the principal value prescription) multiplied with the propagator in the hard 
subprocess (having a causal prescription), 

n eS (v,x) = X ~l. $S(y,s) (37) 
x — y + le 

<f> F (y,x) (38) 



x — y + ie 

$'X(y,x)-7r5(x-y)$ a F (y,x). (39) 



The time- reversal constraint applied to $^(x, y) implies the analogue of Eq. (15), while <& F (x, y) 
has the opposite behavior under time- reversal compared to y). Thus for $^ (x, y) one 

does not have definite behavior under time-reversal symmetry. Specifically, the allowed T- 
even functions of $> F (x,x), G F (x,x) and E F (x,x), can be identified with T-odd functions in 
the effective correlation function $^ efr . This implies that G e A s (x,y) and E A G (x,y) will have an 
imaginary part and this gives rise to two "effective" time-reversal-odd DFs via the imaginary 
part of the e.o.m. 

To say it again in a different way: by partial integration we find for instance 

G A (x,y) = ^-G F (x,y). (40) 
x-y 

If one applies time-reversal invariance, Ga{x,v) will be a real function and Gp{x,y) imaginary. 
So one expects the pole prescription to be the principal value. But when convoluting the pole 
of the matrix element (with the principal value prescription) with the propagator in the hard 
subprocess (with a causal prescription), it is formally possible to shift the imaginary part from the 
pole of the latter to the pole of the former. This will effectively give rise to a causal prescription 
in Eq. fl40"D , instead of a principal value (but without the additional boundary term required 
by time-reversal, cf. Eq. fl3Tp). This implies that G A (x, y) (and also E A (x,y)) will effectively 
acquire an imaginary part. 

For simplicity we neglect intrinsic tranverse momentum, thus we assume ^a(oo)( x ) = ^d(oo)( x )- 

So, by identification we have 



inG F (x,x) = J dylmGf(y,x), (41) 
i7iE F (x,x) = J dylmEf(y,x) (42) 



and then it follows from the e.o.m. that 



xff(x) = i nG F (x,x) = ^ f T(x,S T ), (43) 



xh eS (x) = 2i7rE F {x,x) = ^^Tr[$£(x,x) 7T ^-]. (44) 
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The function receives no gluonic pole contribution, since time-reversal symmetry requires 
H F (x,x) = 0. 

Of course, the mechanism for generating finite projections of <& F (x,x) remains unknown. 
We just can conclude that if there is indeed a non-zero gluonic pole (in the case of non-zero 
antisymmetric boundary conditions), then at twist-three there are two non-zero "effective" T- 
odd DFs, namely fa and h. The first one generates the twist-three single spin asymmetry found 
by Hammon et al. ||, in their notation it is proportional to T(x,x). The second one leads to 
a new asymmetry (see next section). Summarizing, we find for the T-even parametrization of 



ixM 



eff, 



X) It 



ai3 



1 



SrpP + ~h eS (x)^?> 



(45) 



The antisymmetric nonvanishing boundary condition for &a(±oo)( x ) m ight arise from a linear 
A-field, giving a constant field strength (cf. e.g. [17, [H|]). One might also think of an instanton 
background field. In both cases one should interpret infinity to mean 'outside the proton radius'. 
Also, the constant field strength should be understood as an average value of the gluonic chro- 
momagnetic field, which is non-zero due to a correlation with the direction of the proton spin. 
The large distance origin of the asymmetries arising from such a gluonic pole is apparent. 



5 Single spin asymmetries in the Drell-Yan process 

We will now discuss the single spin asymmetries in the Drell-Yan process in case one integrates 
over transverse photon momentum. So one uses the above parametrizations of the correlation 
functions in the expression for the integrated hadron tensor, which after contraction with the 
lepton tensor yields the cross-section. 




Figure 4: Kinematics of the Drell-Yan process in the lepton center of mass frame 



Under the assumption that 



A(oo) 



D(oo) 



we find the following single spin asymmetry 



(hadron- two unpolarized) , given in the lepton center of mass frame: 
4sin(2#) sin(0 Sl ) \S, 



A, 



IT I 



1 + cos 2 9 



Q 



X 



M l x fe(x)f?(x) + M 2 h a x {x)x h d (x) 



£«SA a (*)/i(*). 



(46) 



8 



where 4>s 1 is the angle between S lT and the perpendicular part of the lepton momentum Z, 
lj_ = (g^ — i^i^ + z^z u ^ l v . The first term in the asymmetry (proportional to fx) is the one 
discussed in || (in their notation it is proportional to T(x,x)q(y)), which will also be present in 
DIS (fi(x) = 5(1 — x)). The second term is another, new single spin asymmetry arising in DY 
from a gluonic pole. It is not proportional to T(x, St), but to a chiral-odd projection of $^ in 
the point x = y, cf. Eq. (pPf) . 

6 Conclusions 

We have demonstrated for the Drell-Yan process that the effects of so-called gluonic poles in 
twist-three hadronic matrix elements cannot be distinguished from those of T-odd distribution 
functions. Imaginary phases arising from hard subprocesses together with gluonic poles give 
rise to effective T-odd distribution functions. This leads to single spin asymmetries for the 
Drell-Yan process, such as the one found recently by Hammon et al. ||. We have found a similar 
asymmetry arising from a gluonic pole, which involves chiral odd distribution functions. We have 
moreover shown that the presence of gluonic poles is in accordance with time-reversal invariance 
and requires large distance gluonic fields with antisymmetric boundary conditions. 

We thank A. Schafer for useful discussions. This work was in part supported by the Founda- 
tion for Fundamental Research on Matter (FOM) and the National Organization for Scientific 
Research (NWO). It is also performed in the framework of Grant 96-02-17631 of the Russian 
Foundation for Fundamental Research and Grant iV° 93-1180 from INTAS. 
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